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High transverse momentum (Pt) QCD scattering processes are regarded as a valuable 
tool to study the medium produced in heavy-ion collisions, as due to uncertainty argu- 
ments their cross section should be calculable independent of medium properties whereas 
the medium then modifies only the final state partons emerging from a hard vertex. With 
the heavy-ion physics program at the CERN LHC imminent, the attention of high 
physics in heavy ion collisions is shifting from the observation of hard single hadrons 
to fully reconstructed jets. However, the presence of a background medium at low Pt 
complicates jet-finding as compared to p-p collisions. Monte-Carlo (MC) codes designed 
to simulate the evolution of parton showers evolving into hadron jets are valuable tools 
to understand the complicated interplay between the medium modification of the jet 
and the bias introduced by a specific jet-finding scheme. However, such codes also use 
a set of approximations which needs to be tested against the better understood single 
high Py hadron observables. In this paper, I review the ideas underlying the MC code 
YaJEM (Yet another Jet Energy-loss Model) and present some of the results obtained 
with the code. 



1. The concept of jets in medium 

Jets in vacuum are best understood as a tool to bridge the gap between theory and 
experiment: While using perturbative Quantum Chromodynamics (pQCD) one can 
calculate results in terms of high px parton production cross sections, experimen- 
tally collimated showers of hadrons are observed. In nature, the link between these 
is given by the perturbative parton shower evolution which turns initial high virtu- 
ally partons by braching processes into a shower of softer, lower virtuality partons, 
and by the non-perturbative hadronization of this parton shower. Jets are then 
an operational definition how to treat the observed hadron distribution to 'undo' 
this evolution and combine properties of measured hadrons into an object, the jet, 
which can be compared with calculations on the parton level. For instance, sequen- 
tial recombination algorithms like anti-fc^ try to undo the branchings of the shower 
evolution by combining pairs of hadrons close in phase space. This means that a 
jet is only a meaningful concept in the context of a particular jet finding algorithm, 
and that the choice of the algorithm always corresponds to a bias to pick out a 
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particular subset of all possible evolutions of the original parton. 

In the presence of the soft medium created in heavy-ion collisions, the picture 
becomes more complicated. First, even if a jet which is completely uncorrelated 
with the medium is embedded into the background, the results of jet finding may 
change. This is because the medium constitutes a noisy environment in which the 
soft hadrons of the jet cannot be uniquely identified, thus a medium hadron may be 
combined into part of the jet or a jet hadron may be regarded as part of the back- 
ground. Such complications can be studied by embedding MC jets into a medium 
background and trying to recover them. 

However, in general the role of the medium is much more pronounced. There is 
expected to be significant interaction between the parton shower and the medium, 
i.e. there can be substantial redistribution of energy and momentum from the per- 
turbatively calculable part of the system to the non-perturbative background and 
vice versa: The medium can 'absorb' partons which become sufficiently soft, whereas 
hard scatterings between shower and medium can kick a parton from the medium 
and correlate it with the jet. This has two important consequences: First, the con- 
nection of any jet definition at 'hadron level' with the underlying parton becomes 
tenuous, as a medium hadron which did not originate in the branching process of 
the initial parton may still carry part of its energy and momentum picked up via 
elastic collisions. Instead, the properties of the initial parton are only manifest in the 
flow of energy and momentum carried by both perturbative and non-perturbative 
sector. Second, there is good evidence that the redistribution of energy and mo- 
mentum in the soft bulk is quite different from the dynamics of a parton shower — 
hydrodynamical phenomena such as Shockwaves and diffusion wakes turn out to be 
important Consequently, vacuum jet definitions designed to unfold perturbative 
phenomena may largely miss this part of the evolution. 

There are then two possible classes of jet definitions suitable for the heavy-ion 
environment. The first class focuses on the perturbative part of the evolution and 
eliminates the non-perturbative physics by cuts e.g. on Pp. The medium effect is 
then only apparent from the functional behaviour of jet Raa as a function of the 
cuts as suggested e.g. in El. The second class adapts to the non-perturbative physics 
and tries to capture as much as possible of the initial parton energy and momentum. 
This may require jet definitions very different from the vacuum case. 

With this in mind, some care has to be taken to interpret the results of in- 
medium shower evolution codes correctly. Such codes focus on the perturbative 
side of the evolution only. This means that the code will either make the somewhat 
artificial assumption that no energy and momentum are exchanged between medium 
and parton (as e.g. in Q-PYTHIA^) or that the energy in the parton shower will not 
balance the energy of the initial hard parton because the medium energy balance 
is not explicitly included. 
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2. A description of YaJEM 

YaJEM is an example for an in-medium shower code which permits the exchange 
of energy and momentum between medium and shower , but does not explicitly 
model the medium. The model is described in detail in ISlfilil ft j s based on the 
PYSHOW algorithm IS! which models the shower as a series of 1 — > 2 branchings of 
partons a—* be and to which it reduces in the absence of a medium. The equations 
resulting from pQCD expressions for the branching probabilities are solved by a MC 
in momentum space. In a medium, the medium evolution in position space must be 
connected with the shower evolution equations in momentum space. Within YaJEM 
the link to the medium spacetime dynamics is made by modelling the average time 
for a parton b to branch from parent a given the parton energies and virtualities 
based on the uncertainty relation as 



, V ^b ^b 

whereas the actual time in given branching is generated from the exponential 
branching probability distribution 



(2) 

For simplicity, all partons of a shower are propagated with this time information 
along an eikonal trajectory determined by the shower initiator. This amounts to 
neglecting the spread in transverse space when probing the medium. Currently, 
YaJEM models three different scenarios for the parton-medium interaction, two of 
which modify the kinematics of the propagating parton whereas the last modifies the 
branching probabilities at each vertex (thus, in the first two scenarios the energy of 
the shower is not the energy of the shower initiating parton as there is explicit energy 
transfer between shower and medium, whereas in the last scenario the energy in the 
shower is conserved). For instance, in the RAD scenario, the medium is assumed 
to cause an increase AQq in the virtuality of a parton a based on a local transport 
coefficient as given by the line integral 



P(n) = exp 



n 
(n) 



AQ 2 a = I "dCqiO- (3) 

This modification causes medium-induced radiation. In the DRAG scenario the 
medium is assumed to cause energy and momentum loss along the parton trajec- 
tories, whereas in the FMED scenario the kinematics remains unmodified, but the 
branching probabilities are modified as compared to the vacuum case. In general, 
in the RAD scenario the energy in the final perturbative state is formally larger 
than initially (the depletion of medium partons becoming corelated with the jet is 
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neglected) whereas in the DRAG scenario the final energy is smaller (the absorp- 
tion of partons by the medium us neglected). These effects in principle have to be 
compensated for when embedding YaJEM jets into a background. 

3. Key YaJEM results 

Before applying any in-medium shower evolution to jets, it is useful to understand 
how it relates to known models and observables of leading parton energy loss. An 
important question is to what degree the LPM suppression of subsequent induced 
radiation processes known to lead to the characteristic L 2 dependence of energy 
loss in a constant medium is preserved in a probabilistic MC description. This has 
been computed for YaJEM in ^ and is shown in Fig. Q] left for a 100 GeV charm 
quark propagating through a constant medium in the RAD scenario. 



charm quark, E = 100 GeV, q = 5 GeVVfm 




Fig. 1. Left panel: Mean energy loss per unit pathlcngth from a 100 GeV charm quark propagat- 
ing through a constant medium for two different assumptions about the spacetime picture of the 
shower, compared with an L 2 dependence. Right panel: Examples of medium-modified fragmen- 
tation functions obtained in YaJEM for three different scenarios of parton-medium interaction. 



If only Eq. (fTJ) is used to model the spacetime picture of the shower based on the 
average lifetime of a virtual state, then there is a region in which L 2 dependence 
is seen before finite energy corrections become manifest. However, if Eq. ^ is 
used to determine the lifetime, this region is substantially reduced. Thus, while 
the MC code is capable of preserving this important feature known from analytic 
calculations, an attempt to model more realistically much weakens it. This has also 
been observed in a different framework ®. The disagreement with the pathlength 
dependence deduced from data is currently an unsolved theoretical issue. 

The right panel of Fig. [1] shows the medium- modified fragmentation function in 
all three different scenarios of parton-medium interaction. All lead to a compara- 
ble depletion of the high z region of the fragmentation function, corresponding to 
leading parton energy loss. To see the modification in the low z part more clearly, 
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it is useful to transform into £ = ln(l/z). This is shown in Fig. [21 left panel. 




\ 

Fig. 2. Left panel: The 'hump-backed plateau' distribution dN/di^ obtained in YaJEM for three 
different scenarios of parton-medium interaction. Right panel: The away side per-trigger yield in 
Au-Au collisions normalized to the result in p-p collisions Iaa in 7-h correlations obtained in the 
YaJEM RAD scenario and in a standard leading parton energy loss picture . 



While the DRAG scenario (in which energy flows into the medium) leads to a 
small suppression in this region, induced radiation is manifest as an enhancement of 
the plateau. Experimentally, this region in £ can be probed by 7-hadron correlations. 
The expectation for the away side per trigger yield suppression ratio Iaa from 
YaJEM (with enhanced low z hadron production due to induced radiation) and a 
standard energy loss calculation (where lost energy is assumed to be redistributed 
throughout the whole medium) is shown in Fig. [2] right panel H^l Currently, there 
is no evidence in the data for a rise of Iaa above unity in the data. This might 
be some evidence that a picture of perturbative redistribution of lost energy inside 
the shower only is not justified and that energy flow in the medium needs to be 
accounted for. 

To provide an example for genuine jet observables, Fig. [3] shows the distribution 
of thrust T = max nT ^j P |p" T 1 thrust major T ma j = max nr . n= o ^ ^ ~p an d 

thrust minor T m i n — ^^'^""^ for 100 GeV quark jets propagated through medium 
densities characteristic for LHC Pb-Pb collisions. 

There is a general trend that medium-induced radiation makes the event more 
spherical, regardless if explicit exchange of energy and momentum with the medium 
is modelled or not, but the Pt dependence of the distribution is very different in 
both cases. Thus, assuming that the parton shower has zero momentum exchange 
with the medium may be a bad approximation. This in turn means that the prob- 
lem of energy redistribution by non-perturbative medium degrees of freedom in all 
likelihood will have to be dealt with in order to understand medium-modihed jets. 
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Fig. 3. Distributions of thrust, thrust major and thrust minor for 100 GeV quark jets in vacuum 
and for a medium expected in central 5.5 ATeV Pb-Pb collisions at LHC as computed in YaJEM, 
for all Pt (left) and with a cut of 4 GeV (right). 



4. Conclusions 

Currently medium-modified shower codes have to be regarded as work in progress. 
Qualitatively, some of their results look promising. However, when looking into the 
details, there are often problems in reproducing more differential leading hadron 
observables. Quantitatively, YaJEM awaits comparison with jet measurements at 
RHIC, which, due to the complicated jet finding in a heavy-ion background, is not 
an easy task. There are also some indications that non-perturbative dynamics in 
the medium might play a role for redistributing the energy and momentum of the 
initial hard process in addition to the perturbative dynamics of a parton shower. 
The kinematic range available at the LHC will help much to resolve these questions. 
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